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Abstract: In the article, to bring the system of linear algebraic equations, it is necessary to 
choose such a numbering order that it forms a group of non-zero elements of the matrix. 


The use of the Monte Carlo method leads to a system of linear algebraic equations. To solve the 
system of equations, the modified square root method is used for the symmetric ribbon structure of the 
coefficient matrix since each finite element is connected with a finite number of other finite elements; 
the matrix of the system of equations is always sparsely filled. The location of the coefficients in the 
matrix is closely related to the numbering order of the nodes in the body. 


With random numbering, the non-zero components in the matrix are also randomly distributed, and 
this, in turn, leads to an increase in the calculation time of the problem. 


Therefore, it is necessary to choose such a numbering order, in which non-zero elements are grouped 
near the main diagonal of the matrix, that is, they form a line (ribbon). 


The order of node numbers in the discrete model allows to minimize the width of the tape of non-zero 
coefficients of the equation solving system. Since the coefficients of the matrix of such a system of 
equations and the coefficients of the stiffness matrix of finite elements are symmetric, at the stage of 
solving the system of equations only it is advisable to use diagonal elements and elements. The lower 
triangular matrix is calculated below it. 


In this case, the coefficients of the last matrix also have a band structure. Since the transformation of 
the square root method mainly uses the matrix-to-vector operation, it is necessary to develop a matrix- 
to-vector algorithm for the case where only the coefficients of the lower row of the triangular matrix 
are used. 


if these coefficients are arranged in a row, a rectangular matrix with dimensions Sj, is formed, where n 
is the order of the system of equations, | is non-zero coefficients If these coefficients are arranged in a 
row, a rectangular matrix with dimensions Sj; is formed, where n is the order of the system of 
equations, | is non-zero coefficients. 


Including half the width of the tape. diagonal elements. Moreover, the diagonal elements of the 
original matrix are located in the last Ith column of the Sij matrix. 


To illustrate the changes, let's assume n=9, 1=4. In this case, the lower triangular and rectangular 
matrices Sij have the following forms. 


Published under an exclusive license by open access journals under Volume: 3 Issue: 3 in Mar-2023 
Copyright (c) 2023 Author (s). This is an open-access article distributed under the terms of Creative Commons Attribution 
License (CC BY).To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/ 


128 


8 ae oe ee eT eee 
| J | AET International Journal of Innovative 


Analyses and Emerging Technology 
| e-ISSN: 2792-4025 | http://openaccessjournals.eu | Volume: 3 Issue: 3 


bs 0 0 Os, | 
ty, bop O DO Sy, So, 
ts, by, 33 O $5. 533 S34 
Pay Cay Cag Cag Sa, Sqr S43 S44 
T=|0 ty ty ty ty S=] 85, S52 S53 S54 
0 O tes tea tes bee Soi S62 S63 Seg 
0 OO DO ty ts tye ty Sz, S72 S730 Sag 
0 0 DO DO tes tye ter tes Sg, Sgo S30 Sigg 
0 0 D0 DO DO ty to, tog top Sor So S93 Soa | 


System of coefficients of equations. 


To form the process of multiplaying the matrix sij by the vector xj, the following relation is developed 
and applied, based on the modified coefficients located diagonally and below the matrix: 
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The above relationship makes it possible to perform in the algorithm of the method of square roots 
without zero coefficients and without coefficients of the upper triangular matrix outside the zone of 
non-zero elements using the corresponding diagonal coefficients and coefficients of the lower 
triangular matrix. 
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